Topographical and cell type-specific connectivity of rostral and caudal forelimb corticospinal neuron populations.

Corticospinal neurons (CSNs) synapse directly on spinal neurons, a diverse assortment of cells with unique structural and functional properties necessary for body movements. CSNs modulating forelimb behavior fractionate into caudal forelimb area (CFA) and rostral forelimb area (RFA) motor cortical populations. Despite their prominence, the full diversity of spinal neurons targeted by CFA and RFA CSNs is uncharted. Here, we use anatomical and RNA sequencing methods to show that CSNs synapse onto a remarkably selective group of spinal cell types, favoring inhibitory populations that regulate motoneuron activity and gate sensory feedback. CFA and RFA CSNs target similar spinal neuron types, with notable exceptions that suggest that these populations differ in how they influence behavior. Finally, axon collaterals of CFA and RFA CSNs target similar brain regions yet receive highly divergent inputs. These results detail the rules of CSN connectivity throughout the brain and spinal cord for two regions critical for forelimb behavior.

A patterned human neural tube model using microfluidic gradients.

The human nervous system is a highly complex but organized organ. The foundation of its complexity and organization is laid down during regional patterning of the neural tube, the embryonic precursor to the human nervous system. Historically, studies of neural tube patterning have relied on animal models to uncover underlying principles. Recently, models of neurodevelopment based on human pluripotent stem cells, including neural organoids1-5 and bioengineered neural tube development models6-10, have emerged. However, such models fail to recapitulate neural patterning along both rostral-caudal and dorsal-ventral axes in a three-dimensional tubular geometry, a hallmark of neural tube development. Here we report a human pluripotent stem cell-based, microfluidic neural tube-like structure, the development of which recapitulates several crucial aspects of neural patterning in brain and spinal cord regions and along rostral-caudal and dorsal-ventral axes. This structure was utilized for studying neuronal lineage development, which revealed pre-patterning of axial identities of neural crest progenitors and functional roles of neuromesodermal progenitors and the caudal gene CDX2 in spinal cord and trunk neural crest development. We further developed dorsal-ventral patterned microfluidic forebrain-like structures with spatially segregated dorsal and ventral regions and layered apicobasal cellular organizations that mimic development of the human forebrain pallium and subpallium, respectively. Together, these microfluidics-based neurodevelopment models provide three-dimensional lumenal tissue architectures with in vivo-like spatiotemporal cell differentiation and organization, which will facilitate the study of human neurodevelopment and disease.

A transcriptomic taxonomy of mouse brain-wide spinal projecting neurons.

The brain controls nearly all bodily functions via spinal projecting neurons (SPNs) that carry command signals from the brain to the spinal cord. However, a comprehensive molecular characterization of brain-wide SPNs is still lacking. Here we transcriptionally profiled a total of 65,002 SPNs, identified 76 region-specific SPN types, and mapped these types into a companion atlas of the whole mouse brain1. This taxonomy reveals a three-component organization of SPNs: (1) molecularly homogeneous excitatory SPNs from the cortex, red nucleus and cerebellum with somatotopic spinal terminations suitable for point-to-point communication; (2) heterogeneous populations in the reticular formation with broad spinal termination patterns, suitable for relaying commands related to the activities of the entire spinal cord; and (3) modulatory neurons expressing slow-acting neurotransmitters and/or neuropeptides in the hypothalamus, midbrain and reticular formation for 'gain setting' of brain-spinal signals. In addition, this atlas revealed a LIM homeobox transcription factor code that parcellates the reticulospinal neurons into five molecularly distinct and spatially segregated populations. Finally, we found transcriptional signatures of a subset of SPNs with large soma size and correlated these with fast-firing electrophysiological properties. Together, this study establishes a comprehensive taxonomy of brain-wide SPNs and provides insight into the functional organization of SPNs in mediating brain control of bodily functions.

Spatial atlas of the mouse central nervous system at molecular resolution.

Spatially charting molecular cell types at single-cell resolution across the 3D volume is critical for illustrating the molecular basis of brain anatomy and functions. Single-cell RNA sequencing has profiled molecular cell types in the mouse brain1,2, but cannot capture their spatial organization. Here we used an in situ sequencing method, STARmap PLUS3,4, to profile 1,022 genes in 3D at a voxel size of 194 × 194 × 345 nm3, mapping 1.09 million high-quality cells across the adult mouse brain and spinal cord. We developed computational pipelines to segment, cluster and annotate 230 molecular cell types by single-cell gene expression and 106 molecular tissue regions by spatial niche gene expression. Joint analysis of molecular cell types and molecular tissue regions enabled a systematic molecular spatial cell-type nomenclature and identification of tissue architectures that were undefined in established brain anatomy. To create a transcriptome-wide spatial atlas, we integrated STARmap PLUS measurements with a published single-cell RNA-sequencing atlas1, imputing single-cell expression profiles of 11,844 genes. Finally, we delineated viral tropisms of a brain-wide transgene delivery tool, AAV-PHP.eB5,6. Together, this annotated dataset provides a single-cell resource that integrates the molecular spatial atlas, brain anatomy and the accessibility to genetic manipulation of the mammalian central nervous system.

The neurons that restore walking after paralysis.

A spinal cord injury interrupts pathways from the brain and brainstem that project to the lumbar spinal cord, leading to paralysis. Here we show that spatiotemporal epidural electrical stimulation (EES) of the lumbar spinal cord1-3 applied during neurorehabilitation4,5 (EESREHAB) restored walking in nine individuals with chronic spinal cord injury. This recovery involved a reduction in neuronal activity in the lumbar spinal cord of humans during walking. We hypothesized that this unexpected reduction reflects activity-dependent selection of specific neuronal subpopulations that become essential for a patient to walk after spinal cord injury. To identify these putative neurons, we modelled the technological and therapeutic features underlying EESREHAB in mice. We applied single-nucleus RNA sequencing6-9 and spatial transcriptomics10,11 to the spinal cords of these mice to chart a spatially resolved molecular atlas of recovery from paralysis. We then employed cell type12,13 and spatial prioritization to identify the neurons involved in the recovery of walking. A single population of excitatory interneurons nested within intermediate laminae emerged. Although these neurons are not required for walking before spinal cord injury, we demonstrate that they are essential for the recovery of walking with EES following spinal cord injury. Augmenting the activity of these neurons phenocopied the recovery of walking enabled by EESREHAB, whereas ablating them prevented the recovery of walking that occurs spontaneously after moderate spinal cord injury. We thus identified a recovery-organizing neuronal subpopulation that is necessary and sufficient to regain walking after paralysis. Moreover, our methodology establishes a framework for using molecular cartography to identify the neurons that produce complex behaviours.

Movement is governed by rotational neural dynamics in spinal motor networks.

Although the generation of movements is a fundamental function of the nervous system, the underlying neural principles remain unclear. As flexor and extensor muscle activities alternate during rhythmic movements such as walking, it is often assumed that the responsible neural circuitry is similarly exhibiting alternating activity1. Here we present ensemble recordings of neurons in the lumbar spinal cord that indicate that, rather than alternating, the population is performing a low-dimensional 'rotation' in neural space, in which the neural activity is cycling through all phases continuously during the rhythmic behaviour. The radius of rotation correlates with the intended muscle force, and a perturbation of the low-dimensional trajectory can modify the motor behaviour. As existing models of spinal motor control do not offer an adequate explanation of rotation1,2, we propose a theory of neural generation of movements from which this and other unresolved issues, such as speed regulation, force control and multifunctionalism, are readily explained.

Toll-like receptors 2 and 4 differentially regulate the self-renewal and differentiation of spinal cord neural precursor cells.

BACKGROUND: Toll-like receptors (TLRs) represent critical effectors in the host defense response against various pathogens; however, their known function during development has also highlighted a potential role in cell fate determination and neural differentiation. While glial cells and neural precursor cells (NPCs) of the spinal cord express both TLR2 and TLR4, their influence on self-renewal and cell differentiation remains incompletely described. METHODS: TLR2, TLR4 knock-out and the wild type mice were employed for spinal cord tissue analysis and NPCs isolation at early post-natal stage. Sox2, FoxJ1 and Ki67 expression among others served to identify the undifferentiated and proliferative NPCs; GFAP, Olig2 and ß-III-tubulin markers served to identify astrocytes, oligodendrocytes and neurons respectively after NPC spontaneous differentiation. Multiple comparisons were analyzed using one-way ANOVA, with appropriate corrections such as Tukey's post hoc tests used for comparisons. RESULTS: We discovered that the deletion of TLR2 or TLR4 significantly reduced the number of Sox2-expressing NPCs in the neonatal mouse spinal cord. While TLR2-knockout NPCs displayed enhanced self-renewal, increased proliferation and apoptosis, and delayed neural differentiation, the absence of TLR4 promoted the neural differentiation of NPCs without affecting proliferation, producing long projecting neurons. TLR4 knock-out NPCs showed significantly higher expression of Neurogenin1, that would be involved in the activation of this neurogenic program by a ligand and microenvironment-independent mechanism. Interestingly, the absence of both TLR2 and TLR4, which induces also a significant reduction in the expression of TLR1, in NPCs impeded oligodendrocyte precursor cell maturation to a similar degree. CONCLUSIONS: Our data suggest that Toll-like receptors are needed to maintain Sox2 positive neural progenitors in the spinal cord, however possess distinct regulatory roles in mouse neonatal spinal cord NPCs-while TLR2 and TLR4 play a similar role in oligodendrocytic differentiation, they differentially influence neural differentiation.

Transcription factor Zfp276 drives oligodendroglial differentiation and myelination by switching off the progenitor cell program.

In oligodendrocytes of the vertebrate central nervous system a complex network of transcriptional regulators is required to ensure correct and timely myelination of neuronal axons. Here we identify Zfp276, the only mammalian ZAD-domain containing zinc finger protein, as a transcriptional regulator of oligodendrocyte differentiation and central myelination downstream of Sox10. In the central nervous system, Zfp276 is exclusively expressed in mature oligodendrocytes. Oligodendroglial deletion of Zfp276 led to strongly reduced expression of myelin genes in the early postnatal mouse spinal cord. Retroviral overexpression of Zfp276 in cultured oligodendrocyte precursor cells induced precocious expression of maturation markers and myelin genes, further supporting its role in oligodendroglial differentiation. On the molecular level, Zfp276 directly binds to and represses Sox10-dependent gene regulatory regions of immaturity factors and functionally interacts with the transcriptional repressor Zeb2 to enable fast transition of oligodendrocytes to the myelinating stage.

The histone demethylase Kdm6b regulates subtype diversification of mouse spinal motor neurons during development.

How a single neuronal population diversifies into subtypes with distinct synaptic targets is a fundamental topic in neuroscience whose underlying mechanisms are unclear. Here, we show that the histone H3-lysine 27 demethylase Kdm6b regulates the diversification of motor neurons to distinct subtypes innervating different muscle targets during spinal cord development. In mouse embryonic motor neurons, Kdm6b promotes the medial motor column (MMC) and hypaxial motor column (HMC) fates while inhibiting the lateral motor column (LMC) and preganglionic motor column (PGC) identities. Our single-cell RNA-sequencing analyses reveal the heterogeneity of PGC, LMC, and MMC motor neurons. Further, our single-cell RNA-sequencing data, combined with mouse model studies, demonstrates that Kdm6b acquires cell fate specificity together with the transcription factor complex Isl1-Lhx3. Our study provides mechanistic insight into the gene regulatory network regulating neuronal cell-type diversification and defines a regulatory role of Kdm6b in the generation of motor neuron subtypes in the mouse spinal cord.

Regulation of stem cell identity by miR-200a during spinal cord regeneration.

Axolotls are an important model organism for multiple types of regeneration, including functional spinal cord regeneration. Remarkably, axolotls can repair their spinal cord after a small lesion injury and can also regenerate their entire tail following amputation. Several classical signaling pathways that are used during development are reactivated during regeneration, but how this is regulated remains a mystery. We have previously identified miR-200a as a key factor that promotes successful spinal cord regeneration. Here, using RNA-seq analysis, we discovered that the inhibition of miR-200a results in an upregulation of the classical mesodermal marker brachyury in spinal cord cells after injury. However, these cells still express the neural stem cell marker sox2. In vivo cell tracking allowed us to determine that these cells can give rise to cells of both the neural and mesoderm lineage. Additionally, we found that miR-200a can directly regulate brachyury via a seed sequence in the 3'UTR of the gene. Our data indicate that miR-200a represses mesodermal cell fate after a small lesion injury in the spinal cord when only glial cells and neurons need to be replaced.

Bupivacaine reduces GlyT1 expression by potentiating the p-AMPKα/BDNF signalling pathway in spinal astrocytes of rats.

Bupivacaine, a local anaesthetic, is widely applied in the epidural or subarachnoid space to clinically manage acute and chronic pain. However, the underlying mechanisms are complex and unclear. Glycine transporter 1 (GlyT1) in the spinal cord plays a critical role in various pathologic pain conditions. Therefore, we sought to determine whether bupivacaine exerts its analgesic effect by regulating GlyT1 expression and to determine the underlying mechanisms of regulation. Primary astrocytes prepared from the spinal cord of rats were treated with bupivacaine. The protein levels of GlyT1, brain-derived neurotrophic factor (BDNF) and phosphorylated adenosine 5'-monophosphate (AMP)-activated protein kinase α (p-AMPKα) were measured by western blotting or immunofluorescence. In addition, 7,8-dihydroxyflavone (7,8-DHF, BDNF receptor agonist) and AMPK shRNA were applied to verify the relationship between the regulation of GlyT1 by bupivacaine and the p-AMPKα/BDNF signalling pathway. After treatment with bupivacaine, GlyT1 expression was diminished in a concentration-dependent manner, while the expression of BDNF and p-AMPK was increased. Moreover, 7,8-DHF decreased GlyT1 expression, and AMPK knockdown suppressed the upregulation of BDNF expression by bupivacaine. Finally, we concluded that bupivacaine reduced GlyT1 expression in spinal astrocytes by activating the p-AMPKα/BDNF signalling pathway. These results provide a new mechanism for the analgesic effect of intrathecal bupivacaine in the treatment of acute and chronic pain.

Melatonin ameliorates bupivacaine-induced spinal neurotoxicity in rats by suppressing neuronal NLRP3 inflammasome activation.

Bupivacaine is a common local anesthetic that causes neurotoxicity when used at clinical concentrations. Melatonin (MT), is a potent neuroprotective molecule. The study aimed to characterize the neuroprotective effects of MT on spinal neurotoxicity induced by bupivacaine in rats. It showed that bupivacaine, by intrathecal injection, induced spinal injury, and that the protein levels of Nod-like receptor protein 3 (NLRP3), cleaved caspase-1, and the N-terminal region of gasdermin D (GSDMD-N) were significantly increased. NLRP3 was expressed mainly in neurons and microglia. MT treatment ameliorated bupivacaine-induced spinal cord injury in rats by suppressing activation of neuronal NLRP3 inflammasomes.

Developmental Formation of the GABAergic and Glycinergic Networks in the Mouse Spinal Cord.

Gamma-aminobutyric acid (GABA) and glycine act as inhibitory neurotransmitters. Three types of inhibitory neurons and terminals, GABAergic, GABA/glycine coreleasing, and glycinergic, are orchestrated in the spinal cord neural circuits and play critical roles in regulating pain, locomotive movement, and respiratory rhythms. In this study, we first describe GABAergic and glycinergic transmission and inhibitory networks, consisting of three types of terminals in the mature mouse spinal cord. Second, we describe the developmental formation of GABAergic and glycinergic networks, with a specific focus on the differentiation of neurons, formation of synapses, maturation of removal systems, and changes in their action. GABAergic and glycinergic neurons are derived from the same domains of the ventricular zone. Initially, GABAergic neurons are differentiated, and their axons form synapses. Some of these neurons remain GABAergic in lamina I and II. Many GABAergic neurons convert to a coreleasing state. The coreleasing neurons and terminals remain in the dorsal horn, whereas many ultimately become glycinergic in the ventral horn. During the development of terminals and the transformation from radial glia to astrocytes, GABA and glycine receptor subunit compositions markedly change, removal systems mature, and GABAergic and glycinergic action shifts from excitatory to inhibitory.

Structure of Long-Range Direct and Indirect Spinocerebellar Pathways as Well as Local Spinal Circuits Mediating Proprioception.

Proprioception, the sense of limb and body position, generates a map of the body that is essential for proper motor control, yet we know little about precisely how neurons in proprioceptive pathways are wired. Defining the anatomy of secondary neurons in the spinal cord that integrate and relay proprioceptive and potentially cutaneous information from the periphery to the cerebellum is fundamental to understanding how proprioceptive circuits function. Here, we define the unique anatomic trajectories of long-range direct and indirect spinocerebellar pathways as well as local intersegmental spinal circuits using genetic tools in both male and female mice. We find that Clarke's column neurons, a major contributor to the direct spinocerebellar pathway, has mossy fiber terminals that diversify extensively in the cerebellar cortex with axons terminating bilaterally, but with no significant axon collaterals within the spinal cord, medulla, or cerebellar nuclei. By contrast, we find that two of the indirect pathways, the spino-lateral reticular nucleus and spino-olivary pathways, are in part, derived from cervical Atoh1-lineage neurons, whereas thoracolumbar Atoh1-lineage neurons project mostly locally within the spinal cord. Notably, while cervical and thoracolumbar Atoh1-lineage neurons connect locally with motor neurons, no Clarke's column to motor neuron connections were detected. Together, we define anatomic differences between long-range direct, indirect, and local proprioceptive subcircuits that likely mediate different components of proprioceptive-motor behaviors.SIGNIFICANCE STATEMENT We define the anatomy of long-range direct and indirect spinocerebellar pathways as well as local spinal proprioceptive circuits. We observe that mossy fiber axon terminals of Clarke's column neurons diversify proprioceptive information across granule cells in multiple lobules on both ipsilateral and contralateral sides, sending no significant collaterals within the spinal cord, medulla, or cerebellar nuclei. Strikingly, we find that cervical spinal cord Atoh1-lineage neurons form mainly the indirect spino-lateral reticular nucleus and spino-olivary tracts and thoracolumbar Atoh1-lineage neurons project locally within the spinal cord, whereas only a few Atoh1-lineage neurons form a direct spinocerebellar tract.

Single-cell RNA sequencing reveals Nestin+ active neural stem cells outside the central canal after spinal cord injury.

Neural stem cells (NSCs) in the spinal cord hold great potential for repair after spinal cord injury (SCI). The ependyma in the central canal (CC) region has been considered as the NSCs source in the spinal cord. However, the ependyma function as NSCs after SCI is still under debate. We used Nestin as a marker to isolate potential NSCs and their immediate progeny, and characterized the cells before and after SCI by single-cell RNA-sequencing (scRNA-seq). We identified two subgroups of NSCs: the subgroup located within the CC cannot prime to active NSCs after SCI, while the subgroup located outside the CC were activated and exhibited the active NSCs properties after SCI. We demonstrated the comprehensive dynamic transcriptome of NSCs from quiescent to active NSCs after SCI. This study reveals that Nestin+ cells outside CC were NSCs that activated upon SCI and may thus serve as endogenous NSCs for regenerative treatment of SCI in the future.

New Modified Recombinant Botulinum Neurotoxin Type F with Enhanced Potency.

Botulinum neurotoxins (BoNTs) are notorious toxins and powerful agents and can be lethal, causing botulism, but they are also widely used as therapeutics, particularly to treat neuromuscular disorders. As of today, the commercial BoNT treatments available are from native A or B serotypes. Serotype F has shown efficacy in a clinical trial but has scarcely been used, most likely due to its medium duration of effect. Previously, the uniqueness of the light chain of the F7 subtype was identified and reported, showing an extended interaction with its substrates, VAMPs 1, 2 and 3, and a superior catalytic activity compared to other BoNT/F subtypes. In order to more extensively study the properties of this neurotoxin, we engineered a modified F7 chimera, mrBoNT/F7-1, in which all the regions of the neurotoxin were identical to BoNT/F7 except the activation loop, which was the activation loop from BoNT/F1. Use of the activation loop from BoNT/F1 allowed easier post-translational proteolytic activation of the recombinant protein without otherwise affecting its properties. mrBoNT/F7-1 was expressed, purified and then tested in a suite of in vitro and in vivo assays. mrBoNT/F7-1 was active and showed enhanced potency in comparison to both native and recombinant BoNT/F1. Additionally, the safety profile remained comparable to BoNT/F1 despite the increased potency. This new modified recombinant toxin F7 could be further exploited to develop unique therapeutics to address unmet medical needs.

Purification of astrocyte subtype based on a double reporter mice approach for downstream transcription profiling.

Characterizing the molecular signature of a cell subtype leads to a better understanding of cell diversity, as this molecular data can identify new cellular markers and offer insights about cell function. Here, we describe an efficient protocol to separate a subtype of astrocytes, the Olig2-AS, from other glial cells by using a double reporter mouse approach and to determine the transcriptome profile of the Olig2-AS from the postnatal spinal cord using RNA-sequencing analysis. For complete details on the use and execution of this protocol, please refer to Ohayon et al. (2021).

Trpm5 channels encode bistability of spinal motoneurons and ensure motor control of hindlimbs in mice.

Bistable motoneurons of the spinal cord exhibit warmth-activated plateau potential driven by Na+ and triggered by a brief excitation. The thermoregulating molecular mechanisms of bistability and their role in motor functions remain unknown. Here, we identify thermosensitive Na+-permeable Trpm5 channels as the main molecular players for bistability in mouse motoneurons. Pharmacological, genetic or computational inhibition of Trpm5 occlude bistable-related properties (slow afterdepolarization, windup, plateau potentials) and reduce spinal locomotor outputs while central pattern generators for locomotion operate normally. At cellular level, Trpm5 is activated by a ryanodine-mediated Ca2+ release and turned off by Ca2+ reuptake through the sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) pump. Mice in which Trpm5 is genetically silenced in most lumbar motoneurons develop hindlimb paresis and show difficulties in executing high-demanding locomotor tasks. Overall, by encoding bistability in motoneurons, Trpm5 appears indispensable for producing a postural tone in hindlimbs and amplifying the locomotor output.

A shared transcriptional code orchestrates temporal patterning of the central nervous system.

The molecular mechanisms that produce the full array of neuronal subtypes in the vertebrate nervous system are incompletely understood. Here, we provide evidence of a global temporal patterning program comprising sets of transcription factors that stratifies neurons based on the developmental time at which they are generated. This transcriptional code acts throughout the central nervous system, in parallel to spatial patterning, thereby increasing the diversity of neurons generated along the neuraxis. We further demonstrate that this temporal program operates in stem cell-derived neurons and is under the control of the TGFß signaling pathway. Targeted perturbation of components of the temporal program, Nfia and Nfib, reveals their functional requirement for the generation of late-born neuronal subtypes. Together, our results provide evidence for the existence of a previously unappreciated global temporal transcriptional program of neuronal subtype identity and suggest that the integration of spatial and temporal patterning mechanisms diversifies and organizes neuronal subtypes in the vertebrate nervous system.

TNF-α up-regulates Nanog by activating NF-κB pathway to induce primary rat spinal cord astrocytes dedifferentiation.

AIMS: Astrocytes re-acquire stem cell potential upon inflammation, thereby becoming a promising source of cells for regenerative medicine. Nanog is an essential transcription factor to maintain the characteristics of stem cells. We aimed to investigate the role of Nanog in astrocyte dedifferentiation. MAIN METHODS: TNF-α was used to induce the dedifferentiation of primary rat spinal cord astrocytes. The expression of immature markers CD44 and Musashi-1 was detected by qRT-PCR and immunofluorescence. The Nanog gene is knocked down by small interference RNA. Nanog expression was measured by qRT-PCR and western blotting. BAY 11-7082 was used to suppress NF-κB signals in astrocytes. NF-κB signaling was evaluated by Western blotting. KEY FINDINGS: Our results showed that TNF-α promoted the re-expression of CD44 and Musashi-1 in astrocytes. Dedifferentiated astrocytes could be induced to differentiate into oligodendrocyte lineage cells indicating that the astrocytes had pluripotency. In addition, TNF-α treatment activated NF-κB signaling pathway and up-regulated Nanog. Knockdown of Nanog reversed the increase of CD44 and Musashi-1 induced by TNF-α without affecting the activation of NF-κB signaling. Importantly, blocking NF-κB signaling by BAY 11-7082 inhibited the expression of immature markers suggesting that TNF-α induces dedifferentiation of astrocytes through the NF-κB signaling pathway. BAY 11-7082 could also inhibit the expression of Nanog, which indicated that Nanog was regulated by NF-κB signaling pathway. SIGNIFICANCE: These findings indicate that activation of the NF-κB signaling pathway through TNF-α leads to astrocytes dedifferentiation via Nanog. These results expand our understanding of the mechanism of astrocytes dedifferentiation.